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Ultra-thin 3D topological insulators provide a stage to study the surface physics of such materials by 
minimizing the bulk contribution. Further, the experimentally verified snowflake like structure of the 
Fermi surface leads to a hexagonal warping term, and this discourse examines it as a perturbation 
in the presence of a magnetic field. We find that there are corrections to both energy dispersion 
and eigenstates which in turn alter the density of states in the presence of a magnetic field. Both 
the quantum capacitance and the Hall coefficient are evaluated analytically and it is shown here 
that we recover their established forms along with small corrections which preserve the object 
of treating hexagonal warping perturbatively. In our approach, the established Hall conductivity 
expression develops several minute correction terms and thus its behavior remains largely unaffected 
due to warping. The zero-temperature quantum capacitance exhibits SdH oscillations with reduced 
frequencies, with a lowered average capacitance with increased warping of the Fermi surface, while 
maintaining the usual amplitudes. 


PACS numbers: 73.20.-r, 75.70.Tj, 73.43.Cd 


I. INTRODUCTION 

Recent studies on a new state of matter called a topological insulator (TI) have shown potential advancement for 
spintronic applications [1-3] due to an inherent spin-momentum locking described by the Rashba [4] Hamiltonian, 
characterized by a Dirac piece. Materials like Bi 2 Te 3 , Sb 2 Te 3 , and Bi 2 Se 3 have been known to exhibit the properties 
of a TI [5], which involves conduction of electrons on the surface of an insulator. Scanning tunneling microsopy [6] and 
photoemission spectroscopy [7-9] experiments have shown the existence of the single Dirac cone in TI’s, verifying the 
proposed linear energy spectrum. However, as most theoretical insulators are in practice conductors due to doping 
and crystal lattice defects, transport phenomena are difficult to map accurately, and to this end, research is now 
directed towards thin film TI’s which facilitate the exclusion of contributions from the bulk [10-12]. Essentially, for 
some 3D TIs (e.g. Bi 2 Te 3 ), when one increases chemical potential beyond the Dirac point, a continuous change in 
the Fermi surface can be obtained due to helical Dirac Fermions [14, 15] existing on the surface. The shape of the 
corresponding Fermi surface changes from a circle to a snowflake like structure [13]. This phenomenon is well described 
by a hexagonal warping term proposed by Fu [16] in addition to the original Rashba Hamiltonian. It is worth noting 
that this new Hamiltonian is 3-fold invariant while the corresponding energy dispersion is 6-fold invariant due to time 
reversal symmetry [16]. 

More significantly, in the ultra-thin limit the top and bottom states of a TI begin to hybridize [17, 18] leading to 
a cross over from 3D to 2D surface states. This is easily accounted for by adding a small hybridization term to the 
TI Hamiltonian and it can be shown that Schrodinger piece, which is the generator of particle-hole asymmetry [19], 
is negligible in this context. Properties of ultra-thin TI’s such as magneto-optical response [20, 21] and quantum 
spin Hall [22] and anomalous Hall [23] effects, and possible exitonic super-fluidity [24] have been studied theoretically, 
along with spin texture and circular dichroism in the absence of a magnetic field [25] . A literature search reveals that 
most theoretical studies in this this field have been conducted assuming a uniform Fermi surface and so it is useful to 
study the effects of the inhomogeneity induced by the hexagonal warping term on different physical quantities. We 
consider the impact of the warping term on the Hall conductivity and the zero temperature quantum capacitance in 
ultrathin TI’s, and illucidate on the characteristics of these quantities which are different from the non-warping case 
and maybe experimentally verifiable. 

However, the exact solution to the Hamiltonian in its gauge invariant form in the presence of a magnetic field as 
well as warping is likely mathematically intractable as are most an-harmonic oscillators. To our advantage, we treat 
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the hexagonal warping parameter A to be small [25] with the Dirac piece being the dominant contribution up to 
the Fermi surface. As a consequence, in this paper we will be looking at the effects of hexagonal warping on Hall 
conductivity and quantum capacitance in the presence of a magnetic field using time independent perturbation theory. 
The corrections to the energy dispersion and states of the Rashba Hamiltonian with the Zeeman and hybridization 
pieces will be determined treating the hexagonal warping term perturbatively. 

The following is the structure of the paper. In section H we reiterate the solution to the Rashba Hamiltonian along 
with the Zeeman and hybrization terms in the presence of a magnetic field [26], ignoring the Schrodinger term. We 
then consider the hexagonal warping term perturbatively and we derive the first and second order corrections to the 
energy dispersion. We present an analytical calculation for the asymptotic density of states (DOS) in the presence 
of a magnetic field using the technique of self energy. In section IV we determine the correction to the eigenstates of 
the Hamiltonian stated above using first order perturbation theory and use this result to evaluate the Hall coefficient. 
Quantum capacitance is determined analytically in Section V using the DOS obtained in Section HI. Conclusions are 
presented in Section VI. 


II. THE ULTRATHIN TI HAMILTONIAN AND THE WARPING SPECTRUM 


We begin by considering a Hamiltonian of a TI, similar to the one proposed by Fu, but with an additional piece to 
account for the hybridization of the top and bottom states in the ultra-thin regime [25]. 


H = hVFicTxTTy - T^CTyTT^) + -f AfT^) + ^{11% + TT^) 


( 1 ) 


The first term in the r.h.s of (1) is due to spin-momentum locking effect, with vp as the Fermi velocity, cr^ 
and Ti with i = x,y, z are thel spin and the surface pseudospin degrees of freedom. Here n denotes the minimal 
momentum and = ±1 denotes symmetric and anti-symmetric surface states respectively. The second term has two 
components, the first of which is for Zeeman interaction = sp-bB ^ because of the external magnetic field 

B = Bz with g and gs being the electron g factor and Bohr magneton respectively. The second component Aj is 
the hybridization energy and it accounts for the tunneling between the two surface states. The last term in the right 
hand sight of (I) is the hexagonal warping term, with A as the warping parameter. We denote H = Hq + Hp where 
Hq = hvF{<7^'!Ty — TzUyTTx) + A^cT^ + AfT^ and Hp = ^CT^(7ri]. -I- n^) for convenience. We ignore the warping term for 
the moment and proceed to state the solution to Hq. In order to diagonalize the Hamiltonian Hq, Landau operators 
are introduced as a = + *’’’ 1 /); where I = is the magnetic length and tt = —fV + f A, 

is the gauge invariant momentum. In terms of the ladder operators Hq is reexpressed as 


lUJc. 


Hq = ■^T^((T+a - a-O^) + A^cr^ -f AtT^, 


( 2 ) 


where ujc = vpI is the characteristic frequency analogous to the cyclotron frequency of a usual 2DEG. 
The single-particle eigenstates have the following general form [26]: 


\nas) = u'^p^\n 


1) T, t) + u'^p^\n, T, 1) + 


1) t) + u^biIb, B,l). 


( 3 ) 


|n,T(R),t (i)) represents the n-th Landau Level (LL) eigenstate [26] pertaining to the Hamiltonian in Eq.(2), on 
the top (bottom) surface with spin up (down). Here a = 0,1, s = ±, n = 0,... , 00 , and are the corresponding 
complex coefficients, which can be obtained as 


/„as-h,-S7nas-,*(-Ij Jnas+,Jnas- 


( 4 ) 


where 


/nasi 


1 r^ATT^ 

2 V ^nas 


( 5 ) 


The corresponding state energy is: 


Encs = (-l)“y2a;2n+(A, + sAt)2. (6) 

Since the complete Hamiltonian also contains an anharmonic term Hp, we proceed to the use of time independent 
perturbation theory to determine the corresponding energy corrections. Only the dynamical component of the ladder 
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operators will be utilized for this purpose as it can be readily shown that the contribution from the field is negligible. 
The first order correction to energy collapses trivially as 

Kas = {nas\^a^{TT\ + Trl_)\nas) = 0. (7) 

As is standard, we therefore proceed with second order energy correction calculations, and the general form for this 
is as follows. 


^nois 


n a. s ^ncts 


< nas\Hp\n a s > p 
Er,. — E„' 


( 8 ) 


Note that the summations in the above expression are onl^ on the primed variables. A little algebraic manipulation 
on the determinant of the matrix element < nas\Hp\n as >, yields 


< nas\E[p\n a s > 


= -2,T,t)-A„,+3<-w“4^<n,T,;|n +3,T,;) (9) 

+ In +2,S,t) - <n,B,i \n +3, B,i) 

+ - 4, T, f) - <n,T,i \n' - 3, T, i) 

+ ^„'-3<ljt^“4t < < n,B,i \n - 3,B,i), 


where and = a/ (n' - p){n -p - l)(n' -p-2). 

For a non-zero result, the conditions (n|n ) = i5„„' ^ 0 and s = s are mandatory. On explicitly writing out each 
term and substituting in the second order energy correction formula, it is easy to arrive at the following expression. 


p2 ^ 

Ifi 


Az-|-sAt 1 A 2 -1-sAi 1 

-4^(1-- )4i- 


(n — 2)(n — 1) 
Enas — ^ 


4r^(-l)“+“'(l + 


En 


E. 


n—Zot's 


{ti + 2) (n + 1) 
Enas — ^ra+3a's 


En 


v4(-l)“+“'(l + ^A±^)4l + 

-C/r3.rv.<? 


E „ / 

n— 60 L s 

A^ -l- sAt , 1 


E. 


n-\-3cx's 


)i (10) 


-v4?T3(1 - - 4^)' 

^nas -^n+3a's 


The complete energy spectrum E^^^ is therefore 


T^l _ rp I jp'Z 

^ncts ~ ^nas “r J^nas 


( 11 ) 


Fig.l shows the plot between the energy (the unperturbed energy and the second order correction to the energy) 
on the y axis with the Landau Level index n on the x axis. One can see that the energy corrections for different 
values of A are small, which justifies A being treated perturbatively within the limits of the Fermi surface. Here, we 
have assumed the following [26, 27]: B= 10 T, A^ = 5 meV, At = 3 meV, and these parameters remain unaltered 
throughout unless stated otherwise. Although for a given A there may exist an n for which En < E^, this breakdown 
cannot be realised in any physical material as the required occupation of LL’s far exceeds that of known materials [16]. 
Fig.2 depicts that the second order correction to energy due to perturbation increases quadratically with A, i.e. the 
warping parameter. We will utilize the energy correction derived above in the subsequent sections. 


III. DENSITY OF STATES CALCULATION 

This section contains the formulation of the density of states considering the warped Fermi surface in a perpendicular 
magnetic field. DOS calculations are a fundamental aspect of material science which facilitate the calculation of the 
Hall conductivity and the longitudinal conductivity amongst others. There’s a strong correlation between the DOS 
and the zero temperature quantum capacitance which we shall see in Section V. Here we proceed using the technique 
of self energy. 
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FIG. 1: Perturbative energy correction in comparison to 
LL spectrum with A = 0.2 and 0.4 eV nm^, q = 0 and 

s = 1. 



FIG. 2: A plot of second order energy corrections for 
different values of warping parameter in the range [0,1]. 


Considering the self energy of the system as S {E), one can write the following recursive relation [29, 30] 

1 




^ E-El-E-{Ey 

Fq is the impurity induced LL broadening in eqn(12). The DOS is then expressed in terms of self energy as 


D{E) = Im[ 




( 12 ) 


(13) 


Using the perturbative energy correction derived in the previous section, the modihed energy of the system is E^ = 
En + Ey where we’ve dropped the labels a and s. Drawing the conclusion that << En, and introducing it in eqn 
(12), we arrive at 


E-(u) = rg^ 


E-Er,-^-{E) 




Ei 


[U-U„-S-(U)]2’ 


(14) 


on Taylor expansion. The form of this expression is of great advantage as we are familar with the closed form for the 
first term [28] as shown below. 


r2 V ^ = 


2-kTIE 


where no = 




u-u„-E-(u) [hwlY 

\{E — E“(if)}^ — (Az + sAt)2] is the pole. 


jt{TTno), 


(15) 


The above sum incorporates a factor of two due to the symmetric Filing of states for either value of a, and 
we will extend this to the algebra pertaining to the second term. This term shares a repeated pole with the first term 
and we utilise the residue theorem and a standard summation formula in its evaluation. Then, 


E-{E) 


second term — 


AEElcotinn) + 4E^E^ 




—TTcosec^iirn) + cotlnrij—E'^ 
an 


(16) 


Noting that if <C 1, neglecting higher order terms in E and introducing the factor due to a, lands us at 


[hwcY 


{hw,Y 


COiTTUQ. 


(17) 


Ignoring the self energy contribution to the pole, and using the method in ref. [3] we can extract the total DOS by 
summing over each s branch. 


D{E) = ^ 


DoiE) 


S=±l 


1 + 2 ^ exp{—u{ 


U=1 


27rroE f 


U7r(if^ — (Az + sA() 
{hv 


1_ < 


(hwcf 


2 / ’ 


(18) 
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where no = [E"^ - (A^ + sAt)^] . 

We’ve recovered the asymptotic form (n > 1) of the DOS [27, 28] with a minute additional term which manifests due 
to warping. It turns out that the above expression can be further simplified by reexpressing the infinite sum as follows: 


with 


S = Re 


2t:TqE 
^ 2 


ea:p{—+ ijiu} 


u—0 


Summing the infinite geometric series we obtain, 

ea;p{7^} 


S = Re 


, ,JE^-{A. + sA,r], 


1 — ea;p{—7 ^}Cos/3 
1 — 2exp{—'y‘^}CosP + exp{—2^'^} 


ea;p{7^} — exp{i(3} 

The final expression for the DOS in the presence of a magnetic field then reads 

DiE) = DiE,l) + D{E,-l), 

with D{E, s) defined as 

1 - exp{-C-^nCos{nl^^^^%^) 


D{E,s) = 


DoiE) 


1 - + ea:p{-2(^)2} 


1 - 


El, 

{hwcY 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 


DOS{/J^ 



FIG. 3: Density of states plotted against magnetic field 
varying in the range O.IT to 8T 

Fig. 3 shows how the DOS varies from B = O.IT to i? = 8T and one can see that SdH oscillations set in around 
5.5T which may be a consequence of the asymptotic form of the DOS. Also, the function is singular at i? = 0 and so 
it cannot be used to recover the DOS in the absence of a magnetic field. This is typical of quantum systems where 
the removal of a quantizing field by setting it equal to zero does not yield desired results. Since, in this paper we 
are dedicated to results in the presence of the magnetic field, we will ignore the behavior of the DOS in the absence 
of a magnetic field. In the next two sections we will discuss about the role of the above calculated DOS on different 
physical quantities. 


IV. HALL CONDUCTIVITY 


The Hall coefficient of any conducting material is of interest to condensed matter physics [31, 32], and to this 
end, we present the following calculation which reveals an analytic form for the Hall coefficient of an ultra-thin TI. 
To determine this,we apply the Kubo formula for which one requires the perturbatively corrected eigenstates of the 
Hamiltonian presented in section H. The form of the first order eigenstate correction is 


IV', 


corrected 


) = \r) + 


n' 


< -0" [HpIV^’^ > 

En. — E„I 


\r > 


(24) 
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Since the calculation here involving the evaluation of the matrix element < ■i/'"' \Hp\tjj'^ > bears a stark similarity to 
the energy correction, we’ll state the result directly, but we note that the equation above needs to be modified with 
sums over a and s as a starting point and as usual the requirement s = s holds. 


[ A „,a S A ^,0:8* „,Q: S i A o,CK<S* ^ A „,Q: S 

^n-3^nTt^n-3Tt ^n-2^nTi^n-3Ti ^n-3“nBt^n-3Bt ^ri-2^nBi^n-3Bi 

\u m a/ = \iitxa/o 2 ^ - 


^n—3 


\n — 3 a s)c 


A „, q ; s _ A r,,cx s j_ A ^ _ A ,,ol s 

A ^n+3 “nTt ^n+3Tt ^ri+2 ^nTl ^n-3Tl ^n+3 ^nBt ^n+3Bt ^n+2 ^nBl ^n-\-3Bl 

+ 2^^^- 


, 'S'n+3 

a 

— 171 Qi 5)0 “1“ ^n—3 a s|^ 3 Of 5)o“l“ C^n+3 a s|^“l“ ^ s)q, 

where the CnS are defined as 


\n 3 a 5)0 


(25) 


[ A t,Q!S* „,a 8 4 ^.Ck 8 1 4 ^,Ck<s* „,a 8 J n.cks* „,a s 

^n-3^nTt^n-3Tt ^^^-2 “n-3T^ ^n-3 ^nSt “n-3S^ 

•^71-3 a s — 2 K _ K I 

J-'n—3 

ex 

[ A Tj.CtS*^ „.OC 8 _ 4 0,Ck8* „,Q: 8 J_ ^ n.CkS* ^iCk 8 _ 4 0,Ck8* „,Q: 8 

^n+3 «„Tt “n+3Tt ^"+2 “„T^ “ra-3T^ ^"+3 “nSt “n+3Bt ^"+2 “nB^ “n+3B^ 

'~'n+3 a s = -^ 2_^ 


^n+3 


(26) 


To clarify, the subscript zero denotes the un-corrected eigenstates. The coupling of states due to the introduction of 
the magnetic field removes the complete spin polarization of the n = 0 LL [26], and the Dirac Fermions don’t possess 
a zero mode anomaly. The Kubo-formula is stated below and subsequently we insert the expression from eqn (25) to 
determine the Hall coefficient. Note that the primed variables in the summation below are dummy indices and bear 
no relevance to prior calculations. 


7xy 


2tt 


Im[(nQ;s|T^ 


')x 


(7"^ 71 a s ) X (71 a s r CT \nas 


nexs^n' oc' s' 


71F {Enas'} Bp {En'a's' ) 

{Enas En' a' s'} 


Keeping only leading order correction terms in jCmP since they’re small, and accounting for like terms, 


(27) 


T r/ I Z y\ / I l\ I I I l\ z x\ \ 1 ’^Fi.Enas) np {En> q.! s>) 

-— 2^ Im [(nasjor (T^jn a s )o X (n a s joT cr |nas)oJ-- 

, , , , {-^nas ^n'ex's') 

nexs^n Oi s 


^ XV - 


2 ^2 


+ —— jCn-a a splm [(n — 3Q;s|or^cr^|n'a's')o X (n'a's'jor^cr^jn — 3as)o] 

TT • ^ 


noLS^n' ex' s' 


2 ^2 


cjoe' 


TlF{Enas} Bp {En'a's' ) 

{Enas En' a' s'} 

Bp{Enas} Bp{En'a's'} 


\CnP3 a 8plm [(7i + 3a5|oT'^cr^|7i'a'5')o X (7l'a'5'|or"^a^|71 + 3a5)o] ^ 

“ , , {Enas - En'a's'}^ 


nas^n' a's 


( 28 ) 


On simplification, we get 
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= ^ (2n + l)[n_F( ^nas ) ~ ^F{^n-\-las)] + {^z + 


'^^(^n+las) '^_F(^nas) -i 


n—0 a=0,l s=± 


^n+lc 


+ H H l^"-3 “ + l)/nas+/L2 a's'- " (fL2 c's' + fLs- " 7^2 c.'s'-fLs+) 

n as^a's' 

^ 2 a's') 

(^ncKS ^n—2 a's')^ 

+ H H 1^—3 “ + l)/nas-/L4 a's'+ " «sVL4 a's'-fLs+ " 7^4 a's'+fLs-) 

n as^a.'s' 

^ 77.p’(£'7T,Q,s) TZp’(-E/7 t,_4 ck's') 

(^nas ^n—4: a's'^‘^ 

+ H H l^"+3 “ + 1)7 Ls+7"+4 „'«'- - Ss'(7"+4 c's' + fLs- - fl+i c.'s'-fLs+) 

71 Ois^a's' 

TlFiEnots) '^F{^n-\-4 a's') 


{E na.s a'sO^ 

+ ~~^ ^ ^ |C'„_3 Q s| (s S + l)7nas-7„+2 a's'+ ~ (7„+2 a's'-7nas+ ~ 7„+2 a's'+7nas-) 

71 Q;S^Q;'s' 

^ 7Zi7'(£^7T,Q.s) TZiT" (£^7 t,_|_2 Ol's' ) 

(^nas ^n-\-2 a's'^^ 


(29) 


The first term in the eqn. (29), denotes the contribution to the Dirac Fermion Hall conductivity for the unperturbed 
Hamiltonian and has been derived apriori [26]. The remaining four terms are attributed to the warping correction, 
but are completely dominated by the first term even in the absence of Zeeman splitting and the hybridization term. 
The usual quantum phase transition is observed with increasing magnetic field at Ep = 0 when axy jumps from zero 
to ±e^/27r. 


V. QUANTUM CAPACITANCE 


Transport in TI’s have been studied extensively in comparison to its electrostatic properties, but the latter form an 
integral part of a comprehensive understanding of their electronic properties. Quantum capacitance measurements, 
especially in doped TI’s [33, 34], can reveal the nature of the temperature dependent DOS and enhance our 
understanding of their applicability to spintronics [1]. Such measurements have been reported on carbon nanotubes 
and mono- and bilayer graphene systems [35, 36] and hold sway on the future of device physics, and can lead 
to performance enhanced held effect transistors [37]. This section is dedicated to the discussion of the quantum 
capacitance in an ultra-thin TI at zero temperature. While our formalism does not incorporate a temperature 
dependence, it is unique in that it accounts for the effects of hexagonal warping. 


Quantum capacitance (Cq) is readily dehned as 

Cq = = e^Dr, (30) 

where Ep the Fermi energy and Dp is the temperature dependent DOS, and shown below is a standard form for it. 


Dp 



dE 


df{E - Ep) 
dEp 


(31) 


Here, f{E) is the Fermi-Dirac distribution function, and in the limit of zero temperature, Dp^Q = D{Ep), since the 
Fermi distribution approaches a Heaviside Step function. Refering to E^^ (the index a is set to 1) as Ai?„ in order 























to avoid ambiguity, it is easy to write the expression of D{Ep) from equation (18) as 


D{Ef) = Y, 


DoiE) 


s=±l 


1 - 


1 - 2ea;p{-(gg^)2}Cos(7r ) + exp{-2(^g^)2} 


1 - 


< 

(hwc)' 


By appropriate substitution for Ep a complete expression for zero temperature quantum capacitance 
obtained. 


e2 

(32) 
can be 


Quantum Capacitance{//F cm ) 



Ejr(eF) 


FIG. 4: A plot of quantum capacitance with the Fermi en¬ 
ergy depicting the oscillations around the point of charge 
neutrality, which decrease in amplitude far away from 
Ep = 0 . 



FIG. 5: This shows the ratio of quantum capacitances at 
warping parameter values A = 0.5 and A = 0 versus Fermi 
Energy 



FIG. 6: A plot of quantum capacitance with the magnetic 
field for different values of the warping parameter. 



FIG. 7: This shows the variation of quantum capacitance 
with warping parameter and Fermi energy. 


In Fig.4 we can see oscillations in quantum capacitance with consistent [27, 28] but additional amplitude damping 
as a contribution from the hexagonal warping term, which is small due to the perturbative nature of A. This is 
demonstrated in Fig.5 where the amplitude damping at A = 0.5 eV-nm^ is compared to the case when warping is 
absent by plotting the ratio of quantum capacitances at these values of A. The form of the DOS and consequentially 
that of quantum capacitance is asymptotic, and it fails to recover the electron-hole asymmetry at the point of charge- 
neutrality. However, it proves useful in predicting the behavior of the quantum capacitance for larger values of Fermi 
energy and shows a correspondingly receding envelope of these oscillations. 

SdH oscillations are observed in Fig.6 with known amplitude variation [27] for different values of the warping 
parameter, and with reduced average capacitance and frequency damping as the warping contribution rises, the latter 
of which may be due to the asymptotic form of the DOS. Fig.7 shows the variation of Cq with Ep and A, and we see 
the quantum capacitance increasing with Fermi Energy due to charge accumulation. For a given Ep, we notice that 
the warping term increments lead to a dip in quantum capacitance which is consistent with our earlier results. 

The ratio of values of quantum capacitance (Fig.5) are of the order of 10“^ at magnetic field of 5T, and there is 
little change to this ratio on variation of the magnetic field in the range IT and lOT. A measurement of this ratio upto 
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the stated order, along with observations of reduced frequency SdH oscillations, and receding oscillations in quantum 
capacitance as a function of Fermi energy in different materials such as Bi 2 Te^ and Sb 2 Te^ with experimentally 
established warping parameters [16] in field strength ranges of 2T to 8T might provide experimental proof of the 
results stated here. 


VI. CONCLUSION 

In this paper, we have investigated the properties of a warped ultra-thin TI film in the presence of a perpendicular 
magnetic field. In the ultra-thin limit, hybridization between the top and bottom surfaces of the TI is considered, 
along with the usual relativistic Dirac piece and the term due to Zeeman splitting. In addition and significantly, we’ve 
managed to incorporate the effects of the newly proposed hexagonal warping component by treating it perturbatively, 
and have obtained the correction to the LL’s, demonstrating that the first order energy correction vanishes. The 
second order correction is readily obtained and is shown to be miniscale in comparison to the unperturbed LL’s, 
thereby justifying our methods. 

Calculations for the DOS is done analytically and leads to the recovery of an expression found in a recent pub¬ 
lication [28]. This form is coupled with a factor accounting for the warping contribution and provides a complete 
analytical expressions for the DOS for an ultra-tin TI film. 

Furthermore, the Hall conductivity of the warped TI film is derived by determining the first order state corrections 
and then employing the Kubo formula. Once again, we recover established results [26], including but not limited to 
a quantum phase transition, and go on to show that the warping term has little effect on it, making this a useful 
observation for 2D TI transport. 

Lastly, we’ve turned to the calculation of zero-temperature quantum capacitance, using it’s very direct connection 
to the Fermi surface DOS. In the process of our analysis, we found that quantum capacitance persists with it’s show 
of oscillations as a function of Fermi energy but they are further damped due to the presence of the warping term 
and we predict that these oscillations tend to die out for large values oi Ep. SdH oscillations are also observed with 
remarkable similarity to established results [27] but with a reduced average capacitance and damped frequency. 

While the effect of quantum capacitance can be ignored in most conventional devices, it is crucial to the development 
of nano scale technology which employs ultrathin topological insulators. To this end, it is important to note that 
while experiments have found hexagonal warping, there have been no experiments which have measured the effects 
of hexagonal warping on quantum capacitance in ultrathin TI’s. Experiments have measured Hall conductivity and 
quantum capacitance [38, 39], wherein the sample was thick enough to minimize the effects of hybridization of the 
top and bottom surfaces, and without considering the effects of hexagonal warping. We’ve shown that our results 
are consistent with the findings of such experiments and the corresponding theoretical analyses. As such, one can 
measure the quantum capacitance by varying the channel voltage between a gate and an ultrathin TI, and noting the 
corresponding charge response. We propose that future experiments be conducted on gated devices in conjunction 
with ultrathin TI’s, with parametric ranges consistent with those used in this paper, the results of which may be used 
to validate our claims. 

During this discourse, we have stressed on the derivations of quantities fundamental to almost any condensed 
matter system, both in transport and in electrostatics, and stated possible ranges of parameters to facilitate 
experimental validation. It is felt that significant contributions in the understanding of thin-film TI’s lie in the 
understanding of the hexagonal warping and the results presented here may be adopted to analyse other interesting 
properties of ultrathin topological insulators. 
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